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We report dynamical calculations for large-g structure functions of liquid 

o 

O ■ 4 He at T=1.6 and 2.3 K and compare those with recent MARI data. We 



extend those calculations far beyond the experimental range q < 29k m 



order to study the approach of the response to its asymptotic limit for a 
system with interactions having a strong short-range repulsion. We find only 
small deviations from theoretical \ jq behavior, valid for smooth V. We repeat 
an extraction by Glyde et al of cumulant coefficients from data. We argue that 
fits determine the single atom momentum distribution, but express doubt as 
to the extraction of meaningful Final State Interaction parameters. 
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I. INTRODUCTION 



In the following we discuss different aspects of the response of liquid 4 He. First we report 
predictions which are compared with most recent data. Next, we compute the response for 
q ^ 300A 1 in order to study how Final State Interaction (FSI) effects vanish for those large 
q. In the end we present results of a new model-independent cumulant analysis of data in 
order to extract the single-atom momentum distribution and interaction parameters. 

We address below recent precision data for temperatures below and above the transition 
temperature T c , which have been taken at the Rutherford ISIS facility by means of the MARI 
spectrometer. Those by Andersen etal span neutron momentum transfers 3 < q(A 1 ) < 10 
for T=1.42 K and 3 < g(A _1 ) < 17 for T=2.5 kB, while Azuah's measurements covered 
10 < g(A _1 ) < 29 for T =1.6, 2.3§. The above expand in scope previous information taken 
a few years ago at the IPNS facility at Argonne for q < 23.1 A 1 i. 

From inclusive cross sections one extracts the linear response. The latter is a function of 
two parameters q and u, which in the scattering experiment are the momentum and energy, 
transferred from the projectile to the target. For medium and large q those responses 
contain information on the target, such as the momentum distribution of the constituents 
and prescribed manifestations of their interaction. 

A dynamical calculation of the response requires as input the atom-atom interaction 
and groundstate information, which for the above g-regime are primarily the single-atom 
momentum distribution n(p) and the semi-diagonal two-particle density matrix. 

Using variations of much the same theory, calculations have been made before for 
medium-g, as well as for the higher-g Argonne dataihi. To our knowledge no ab initio 
calculations of the MARI data have been performed until now. 

The same MARI data have recently been approached in an entirely different fashion 
with the purpose to determine in a model-independent way the dominant coefficients in the 
cumulant expansions of the asymptotic and FSI parts of the responselE. Good fits to the 
data were obtained, but those have little in common with dynamic calculations. The latter 
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use as input n(p), additional groundstate information and V whereas, ideally, from cumulant 
fits n(p) and properties of V are extracted. 

As a major result of the above analysis, Glyde etal report the reconstruction of the 
single-atom momentum distribution n(p) in good agreement with accurate theoretical 



lant coefficient which, dependent on the analysis, is reported to be less than 0.65 times the 
calculated value. One then wonders whether the apparent partial fit may have consequences 
on the precision of the reconstructed n(p). 

The following program emerges from the above observations. In Section II we outline 
an approach to high-g responses. In Section III we report computations of the high-g mea- 
surements using the MARI spectrometer and compare those predictions with the data. In 
addition we interpret responses computed out to very high q < 300 A \ The results enable 
the study of the approach of the response to its asymptotic limit for systems with a strong 
short-range repulsion in the interaction between the constituents. In Section IV we present 
fits for cumulant parameters for T=2.3K and compare those with similar results by Glyde 
et a We discuss the discrepancy between the calculated and the extracted FSI parame- 
ters and attribute it to the truncation of the cumulant series. In the conclusion we estimate 
that both experimental and theoretical studies of the response of liquid 4 He at high q may 
have reached a degree of sophistication, beyond which there is little prospect to gain new 
information. 

II. DESCRIPTIONS OF THE LINEAR RESPONSE FOR HIGH Q. 

Consider for infinitely extended liquid 4 He the response per atom in the form 




However, a less satisfactory feat 



ure is the extracted dominant FSI cumu- 




<f>(q,y) = {q/M)S(q,u), 



(1) 



with M the mass of a 4 He atom. p q {t) above is the density operator 
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p q (t) = e- im p q (0)e im 

p 9 (o) = E^ (0) (2) 

j 

Strictly speaking, the symbol (0|...|0) should stand for a canonical average at given T, but 
we shall use instead the ground state in conjunction with T-dependent quantities. 

In the last line of (1) we introduce the reduced response <f)(q,y) with the energy loss u, 
replaced by an alternative kinematic variable y = y(q, cu)0lll 

* = (3) 



q V 2M, 

Upon substitution of (2) into (1) one generates two components of the response. In the 
incoherent part, one tracks the same particle when propagating in the medium, while in the 
coherent part one transfers momentum and energy to a particle distinct from the struck one. 
For q > 8A 1 the response is dominated by the incoherent part and the coherent part one 
can be safely disregarded. 

For the description of the large-g response we shall exploit the theory of Gersch, Ro- 
driguez and Smith (GRS) for smooth interactions V which leads to the following expansion 
for the reduced response in inverse powers of q or of the recoil velocity v q = q/ME (We use 
units h = c = 1 causing all quantities to have dimensions of powers of A or A X ) 

oo / 1 \ n 

= E - ( 4a ) 

n=0 Vt V 

roo 

F (y) = Mm (j>(q,y) = (Air 2 )' 1 / dppn(p) (4b) 
i~ *°° J\y\ 

-F 1 (y) = t(2np)- 1 / ds e iys / drp 2 (r, 0; r - sq, 0)x(q; r, s) (4c) 

f ' ds'Vir - sq) - sV(r - sq) 
Jo 



x(q,r,s)-- 



... etc. 



(4d) 



The function x i n O) resembles an eikonal phase. It differs from it because the integration 
limits on the line integral over the first component of V are not (— oo, oo), as is appropriate 
for on-shell scattering. The finite limits correspond to off-shell scattering described in the 
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coordinate representation. Moreover a second interaction is implicit in (fh|). In the following 
we shall allude to the total expression (|4d|) as the generalized eikonal phase. 

We recall the interpretation of the lowest order terms. For sufficiently large momentum 
transfer q, an atom with initial momentum p recoils with p' = \p + q\ « q ^> (p 2 ) 1 ^ 2 , which 
is larger than the average momentum of an atom in the medium and is moreover in excess of 
any inverse length in the system. The recoiling atom moves therefore too fast to be affected 
by atom-atom collisions and the response is the asymptotic limit Fo(y) for q,u — > oo at 
fixed y. Eq. (flbf) shows its expression in terms of the single-atom momentum distribution, 
normalized as / dp / (27r) 3 n(p) = 1. 

Although the GRS theory is not a perturbation theory in the interaction V, the second 
term in the series (|4a|), linear in V, is entirely due to binary collisions (BC) between the hit 
and any other atom. It accounts for the dominant FSI collecting all contributions oc 1/q. 
This is achieved at the price of introducing the semi-diagonal two-particle density matrix p 2 



in (He 



In another publication Gersch and co-workers suggested an alternative representation for 
the reduced response^ 

<P(q, y) = J dy'F (y - y')R(q, y') = J ^- n (p)R(q, y - p z ) (5a) 

4>{q t8 )= / dye-^<f>(q,y) = Yl ~ 

= F (s)R(q, s) = F (s)exp[n(q, s)) (5b) 

In Eq. (|5a| ) the response is written as a convolution of its asymptotic limit and a FSI factor 
R(q, y). It is frequently convenient to use Fourier transforms Fi(q, s), R(q, s).... In particular 
for the first two terms in (5b) one has (cf. (fjb])-(f4c|)) 

F ( s) = ^0) = |_*_ e - ip4 ,„ w (6a) 

F i( s ) = ~ drp 2 {r - sq, 0; r, 0)x{q; r, s), (6b) 
P J 

with pi(s, 0) = pi(r — sq,r), the single-atom density matrix and p = px(r,r), the number 
density. 



We shall restrict ourselves below to various descriptions of FSI due to BC, starting from 
the corresponding cumulant form (|515| ) and using (|6a|)El 

0(9, s) = — i? 2 (g, s) = — - — exp[fi 2 (g, s)] 

n 2 (?, s) = ip J dr( 2 (r, s)u 2 (q, r, s), (7) 

where Rzfa, s) and fi 2 (g, s) are the BC approximation to the corresponding quantities defined 
in (|5aD and (^). C2 above is defined by 

<M = P * (r -f °' r ' 0) (8a) 
Ur,0) = g 2 (r) (8b) 

with g 2 the pair-distribution function. 

Eq. (0) is the most general cumulant form in the BC approximation for the FSI phase 
fi 2 (g, s) = ln[i? 2 (g, s)}, and distinguishes through uj 2 between different dynamical approaches. 
For instance for smooth interactions V, which would allow for an expansion of the expo- 
nential in (|7|), comparison of ([?]), ( [5b]) and (4d) shows uo 2 to be the generalized eikonal 



phase 

U2,v{q, r, s) = x(q, r, s) (9) 

For interactions with a strong short-range repulsion, the line integral over V in the (off- 
shell) phase (icl) which enters the dominant BC FSI contribution i*i(s), Eq. (|6b|), may 



produce large and even divergent integrals. The standard method to tackle those difficulties 
is by partial summation of selected higher order terms 

iu 2 y -> iu)2,t = e lx ~ 1, (10) 

which amounts to replacing the 'bare' V by a g-dependent effective interaction V — > 
V e ff(q) = t(q), the latter being the off-shell t-matrix, in turn generated by V. Moreover the 
propagation in between collisions is described in the eikonal approximation!!!!. 
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In an alternative regularization for an atom-atom interaction with a strong short- 
range repulsion, one replaces the generalized eikonal phase (Bel) by a semi-classical 



approximation! 1 ^ 1 



aj 2 ,sc(q,r,s) =iq [ ds' 
Jo 



/ 2m Tr , / 2m Tr , 



(11) 



For (2m/q 2 V) <C 1, u> 2> sc coincides with u 2 y, Eq. @. However, in classically forbidden 
regions {2m/q 2 )V > 1, uj 2 ^ sc describes damping, as the dominant imaginary part of V e ff(q) 
in fllOD is expected to do. This will be borne out by calculations. 

Whereas u 2 y is strictly proportional to 1/q, this is no more the case for u 2)t after the 
replacement V — > V e ff(q). The above manifestly introduces g-dependence in coefficients of 
the GRS series (4a), (5b) and in particular in the BC approximation. Taking the latter in 
the cumulant form (|7]) adds to the blurring of the original 1/q dependence. This raises the 
question how the response approaches its asymptotic limit. 

We start with a theoretical analysis of the first cause of additional g-dependence and 
focus on 4 He- 4 He scattering for high lab momenta q. The latter is of distinct diffractive 
nature, typical for interactions with a strong, short-range repulsion. For those, the dominant 
imaginary part of the on-shell scattering amplitude t(q) ~ Im/(g) oc iqcr^ 01 , where the total 
4 He- 4 He cross section a q varies much slower than q itselfi. 

Without entering into details, we state that the off-shell t = V e ff in u 2 , Eq. (0), can 
approximately be related to the on-shell scattering amplitude for elastic scattering. (See 
ref.@ for a more extensive treatment of the parallel discussion for atomic nuclei). It can 
then be shown that the rigorous proportionality of the dominant BC FSI phase Cl 2 y ocl/q 
for a smooth, bare V still holds approximately for f^,*- 

Additional g-dependence is due to the use of the cumulant representation (|7|) but it will 
be small to the extent that FSI are. In conclusion, the reduced response described by (5) 
and (7) is expected to approximately preserve the 1/q signature of the dominant binary 
collision contribution. We shall return below to a numerical confirmation. 
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III. DYNAMICAL CALCULATIONS OF SELECTED MARI 4 HE DATA 



We first mention and discuss the input elements which suffice for the BC approximation 
in any of the forms described in Section II. 

a) The atom-atom interaction Vaziz^- 

b) The single-atom momentum distribution n(p,T) 



n(p; T) = (2iry'S(p)n (T) + (1 - no(T))n""(p; T) 
e^ = M T)Hl-MT))£&?®- (12) 

n o(T < T c ) is the fraction of atoms in the condensed state!!, n NO (p;T) and pi(s,0;T)/p 
above are respectively, the momentum distribution of the normal (uncondensed) atoms and 
its Fourier transform. Path integral Monte Carlo (PIMC) calculations have shown moderate 
T-dependence of n NO (p; T) for T < 4K@0. 

c) The least accessible ground state property required in the BC approximation is the 
semi-diagonal, two-body density matrix which weights the dominant BC FSI terms in fljc| ) 
or ([5b]). Calculations based on a variationally determined groundstate wave function in the 
Hyper- Netted Chain (HNC) formalism produce for £ 2 , Eq. 

(? CN {r, s; f) = g wd (r)g wd (\r - sq\)exp [A(r, s)} 
~ 9wd(r)g wd (\r - sq\)exp [£A 4 (r, s)] 
A 4 (r, s)= P J dr'^gUW - sq\) - l] [^,(r') - l] \g{\r' - r\) - l] (13) 

gwd{f) is an auxiliary function related to what in HNC formalism is called a form factoiS. 
The function A(r, s) formally adds all so-called Abe diagrams and is approximated in (13) 
by the 4-body Abe diagram A 4 (r, s), using in addition a scaling parameter £§E1. 
Far less sophisticated and simpler is the GRS approximationB 



(? RS (r,s) = ^g(r)g(r-sq), (14) 
which interpolates £2 between s = and the Hartree limit for large s0. 



Both options have drawbacks and fail for instance the extended unitarity test 

J drp 2 (r - sq, 0; r, 0) = (A - l) Pl (s, 0) (15) 
which can be written as 

S(s)=p|dr(l-C 2 (r;s))=l (16) 

Using a typical pair distribution function g(r), S(s) above computed with Eq. 
([14]) , produces values up to 1.7 for s=2.0 instead of the exact value 1.0, independent of 
s0. In the HNC case, approximations involved in the evaluation of the Abe terms (|13|) are 
responsible for similar deviations of S(s) from 1. The violation of condition (Jl6|) is intrinsic 
in the GRS approximation (0), no matter what g(r) is used. 

Another important constraint is the fact that the diagonal two-body density matrix 
should coincide with the pair distribution function: ( 2 (r,0) = g(r) (|8b|) . While the GRS 
approximation fulfils that condition by construction, a full evaluation of the Abe terms is 
necessary in the HNC formalism. Demanding the boundary value condition to be fulfilled 
in the mean, one determines the so-called scaling £ by minimizing the following quantity 

a(£) = Jdr\C HNC (r,0;O-g(r)\ 2 

A particular choice of (2 presumably matters for medium q, but for increasing q > 20A 
FSI contributions decrease in importance relative to the asymptotic response. A few-percent 
spread, due to uncertainty in the choice of (2, in anyhow small FSI terms will go unnoticed. 
We thus opted for the expression ([14]) which is numerically much easier to handle than ( |13| ) . 

d) Finally, for a comparison of actual data with predictions the latter have to be folded 
into the experimental resolution function (ER) E(q,y;T) of the instrument. The E(q,y;T) 
corresponding to the q > 20A 1 MARI data are given in Azuah's thesisi and have been 
fitted to the sum of two off-center gaussians. No ER, pertinent to lower g's were available 
to us, thus precluding an analysis for q < 20 A . 

Until this point we did not specify the T-dependence of the theoretical responses. In 
fact one ought to employ quantities computed for given T. In fact there exist experimental 
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dataH and also PIMC studies§!Il on the T dependence of the pair-distribution function 
g(r,T). However, in view of the above arguments we shall use the one for T = 0. By 
the same token ( 2 , Eq. (|H]), and consequently FSI effects will be independent of T. This 
leaves the single-particle density matrices, or equivalently the momentum distributions as 
the only T— dependent quantities in the present analysis. We took n (T = 1.6K) = 0.087 
and pf °(0, s;T = 1.6K) = pi(0, s; T = 2.3K) from calculations for T = 1.54 and 2.5K0R 
The expression for the predicted response is therefore 

<f>{q, y;T>T c ) = J j0^ n (P^ T)R(q, y - p z ) 

<f>(q, y-T< T c ) = n Q (T)R(q, y) + [1 - n (T)} J ^n NO (p; T)R(q, y-p z ), (17) 
which in order to enable a comparison with data, has to be folded into ER 

/oo roc _ 

dy' E(q, y - y'; T)<j>{q, y'; T) = (2^)^ / ds e^ s E(q, s, T)0(g, s; T) (18) 
-oo J —oo 

For future reference we emphasize here that the FSI factor R is from (0) and (|]) seen to be 
independent of the single particle density matrix pi(s,0)/p. In particular for all but pure 
hard-core interactions 

Y\mR(q.s) = l. (19) 

q— >oo 

We thus computed the reduced response (f) E (q,y,T), for the q = 21, 23, 25, 29 A sample 
out of the MARI data. In view of the steady decrease of FSI, this g-range and steps seems 
to be sufficient for our study. We emphasize in particular the case q = 23 A , considered 
because it is the largest q in the older Argonne data setM and for it we shall compare our 
results with others. 

We start with a comparison of our predictions for T = 2.3K and the corresponding data! 
(Figs. 1 a-d). The overall agreement is very good. One notices that, whereas the central 
value for the theoretical response hardly changes for 21 < q(A 1 ) < 29, the data for the 
same, folded in the ER, (f) E (q,0) show g-dependence present in E(q,y). 
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The agreement for T = 1.6K (see Figs. 2 a-d) is slightly worse. The slight staggering 
in the central region for q = 2lA 1 is probably of instrumental origin, but contrary to the 
T = 2.3K case, differences in E(q, y) for q = 21, 29A 1 do not explain the small discrepancies 
in their central regions. We recall that exactly the same input is used as for T = 2.3K and 
that the only extra parameter is the condensate fraction n (T = 1.6K). 

We now reach our second topic. In spite of the fact that no data exist for q > 29A , 
we have extended calculations up to q = 300A . The purpose of the exercise is to obtain 
theoretical information on the approach of the response to its asymptotic limit. 

In Fig. 3a we present the part cj) even (q, y,T = 2.3K) of the predicted response, computed 
with (|7|), (TOP which is even in y. Even in the wings out to y ~ 3.5A 1 those coincide within 
1% amongst themselves and with the asymptotic response F (y), Eq. (ff§). The above 
appears hardly changed, when predictions for 20 < q(A ) < 29 are included: only in the 
immediate neighborhood of y = 0, is there a < 2.5% difference. 

In Fig. 3b we show q<ft odd (q, y,T — 2.3K) which contains the part of the response which is 
odd in y. Having included the q signature of the dominant FS some residual q dependence is 
apparent in the extrema as well as in the wings, but the conclusion is clear: Neither the strong 
short-range repulsion in the atom- atom interaction which forces the use of V e ff(q) = t(q), nor 
the effect of the cumulant representation, much changes the 1/q signature of the dominant 
FSI term in the GRS series ( |4aD for smooth V. The above agrees with our arguments in 
Section II and with our previous results!. 

All reported predictions are based on the use of the t-matrix, i.e. on fllPl). In Section II 
we also mentioned a semi-classical approximation (]TTj) for FSI and found that, except for 
small s, there are considerable differences between the BC phases, calculated by means of 
(0) and 

Ultimately excellent agreement is obtained between the corresponding responses, com- 
puted with (||) and (0). Clearly both the t-matrix and the semi-classical method accurately 
describe the binary collision phase in the salient region just inside the classically forbidden 
region. Contributions from deeper penetration distances are strongly suppressed. 
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We conclude this Section, by a discussion and comparison of predictions for q = 23A 
by other authors. Since the various studies refer to different T and data have been taken at 
different instruments, the natural quantity to compare is the FSI factor R(q, y) assumed to 
be T- independent. 

We start with predictions by Mazzanti et all! which are based on exactly the same BC 
approximation in the V e ff version (|T0|), employing however the variationally derived C% , 
Eq. ([131) . Next we mention Silveii who used, what amounts to the cumulant form (|7]) with 
u>2 — > u)2,t and (2 — > #2, the pair distribution function. In his Hard Core Perturbation Theory 
he disregarded the second part of the total phase \i Eq. (f^dj), which is only permissible 
for a pure hard-core interaction. However, Silver actually constructed the off-shell t-matrix 
in (|TU|), corresponding to the first part in (id) from JWKB partial wave phase shifts for 



a realistic V, which in addition to strong short-range repulsion also included attractive 
components. Nevertheless, he neglected the second component in (4) which does not vanish 
when an attraction is present. We conclude with a path- integral method by Carraro and 
Koonin, who computed high-g FSI using a fixed scattering approximation for the entire 
system with a large, finite number of atoms9. The method requires the parallel calculation 
of the ground state wave function in order to construct the iV-body density matrix, diagonal 
except for one particle, with iV the number of atoms in the sample and which averages the 
response for fixed scatterers. 

Results for R(q, y) for all cases discussed are assembled in Fig. 4 and show occasionally 
substantial differences. However, those are considerably smoothened by the single particle 
momentum distribution (Eq. (gD) or density matrix in and ultimately produce quite 
similar responses!!!. 

IV. CUMULANT EXPANSIONS OF THE RESPONSE. 

We consider below a method which has extensively been applied in the pastS before the 
rediscovery of the 1/q GRS expansion of the response ( Recently it has been brought 
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to the fore again in an attempt to parameterize data without the intervention of a theory. 
The method uses cumulant expansions of the Fourier transforms of the separate asymptotic 



and FSI parts of the response 
response^ 



), with coefficient functions related to moments of the 



(-is) , 



n>2 



n>2 



nl 



(j)(q, s) = exp 
F (s) = exp 
R(q, s) = exp[f2(g, s)] = exp 

L n>3 ' t: J 

Using (|5l]) the various coefficient functions are related by Jl n (q) = a n + (3 n (q)^- Data for 
the response are then compared with the parameterization (p0|) 



(20a) 
(20b) 
(20c) 



<!> e E P (q,y) H dse^Ms, [fl n (q)}), (21) 

7T Jo 

where, as in (|18|) the right-hand side in ( j20a| ) includes ER. In principle, no a priori knowl- 
edge, or even meaning of the cumulant coefficient functions [i n {q) is necessary for a search. 
However, it is natural to use motivated initial values such as calculated ones. 

Consider the normal fluid, in which case the g-independent cumulant coefficients a n 
originating in the asymptotic part of 4>(q, s) can simply be expressed in terms of averages of 
even powers of the momentum of an atom, e.g. 

1 



«2 = ^{P 2 ) 



(22) 



The g-dependent coefficients f3 n (q) relate to the FSI factors R(q, s), f2(g, s) in (^b|) and may 
be written as§ 



Pn(q) 



[(n-l)/2] 



m=l 



n—2m 



[(n-l)/2] 



2 m 



( g* 



m=l 



n—2m 



(23) 



t -l 



For convenience we use q*, the momentum transfer parameter q in units of 10 A . Eq. 
displays q— dependence and defines coefficients z nt2m , which have the same dimensions and 
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may be expressed in the same units as /3 n . The above are operators for dynamical variables 
of the system, averaged over diagonal /-body density matrices and their derivatives of order 
/ < n. For the lowest order cumulant coefficient functions one has 

fo(q) = 61im[Imft 2 (g,s)/s 3 ] (24a) 

1 - 2l 



= 7-< V V ) P2 ( 24b ) 

fa(q) = 241im [Reft 2 (g, s)/s 4 } (24c) 
k(q) = -P 5 4+(-) 3 p 5 2 (24d) 

Vg \V q J 

Contrary to the GRS series in 1/q, with coefficients depending on nearly-diagonal n-particle 
density matrices p n {r\ — sq, r 2 , ■ ■■T n ; n, ...r n ), the moment approach underlying the cumu- 
lant expansion does not produce a systematic g-dependence of the coefficient functions. For 
instance, all j3 n (q) with odd n contain FSI contributions oc 1/q. Eqs. ( |24a|) and ( |24d|) 



illustrate this for the dominant FSI coefficient functions j3 n (q). Thus ^(q) oc 1/q draws 
exclusively on the BC contribution (|7|), while the two components of f3§ are proportional 
to 1/q and (1/q) 3 due to binary, respectively higher order collision contributions, etc. It 
is the expansion of the semi-diagonal p 2 (ri — sq, rj; r\, rj) in s which produces an infinite 
number of contributions /5 2n+ i(g), all of which have parts oc 1/q with coefficients depending 
on diagonal densities (cf. ( [24b|) ) and their derivatives. 

We start with the threshold behavior of the computed BC FSI phase Cl 2 (q,s), Eqs. 
(0), ([TOD and in particular of its imaginary part which, according to (|24a|) , produces the 



dominant FSI parameter f3^(q). We checked that, within a few percent glmf2 2 (g, s) is out 
to s ~ 0.8A independent of q. In particular we could extract the theoretical threshold 
value Z3 = q*/3^ 2 (q) = 6q* lim s _^o[I m ^2(o , ) s)/s 3 ] = 0.555A 3 which over the entire range 
q < 300A 1 is to better than 1% independent of q. For the above reason one cannot 
determine the next order odd-ra coefficients z§ with reasonable precision. 

Regarding parameters of even order in n one finds from Re fi 2 (g, s), Eq. fl24cQ for the 
leading coefficient Z42 = q* 2 P^(q) = — 2.26A 4 . It appears that for even n, the next-to- leading 
^64 ~ q* 2 P6i(q) ~ — 31 A 6 is reasonably well defined. 
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We return to Eq. ( |24b| ) which seems to provide an independent way to calculate /3s(q). 
However, one can actually derive it from ( |24a| ) , using (0) for Cl 2 in either version (|9|) or 
(|T0|) for the generalized eikonal phase u(q, s). It holds for arbitrary semi-diagonal p 2 which 
exactly satisfies fl3l5|)ii. 

Using the same g{r) as in (0) we compute q*/3j 2v (q) = 0.56A 3 . The agreement 

between q*j3^ 2y and q*(3^ is very good, especially in view of the sensitivity of q*P^ 2y on the 
precise shape of g(r) where the Laplacian of Vis large. The extraordinary stability of the 
extracted q*f3% confirms the numerical consistency of the calculation. 

As has been mentioned before, all FSI functions have been assumed to be T-independent 
and we have used values for T = 0. In order to estimate the influence of the temperature 
we have also calculated q*P% v using a g{r) obtained with PIMC at T = The result 

q*j3j 2y (T = 2.8K) = 0A7A~ 3 confirms its sensitivity on the precise shape of g(r). 

We now report several results for cumulant fits ( pi]) , obtained with the CERN MINUIT 
code and applied to the ten T = 2.3K data sets in the range 21 < q(A ) < 29 from 
Azuah's thesisS We first note that the integrated strengths of the data there appear for all 
q to be approximately 1.4% in excess of the exact result 1. By construction that demand 
is exactly fulfilled by a cumulant expansion fl20ap in, no matter what approximation. We 



also considered data cut at y « 3.0 — 3.4A 1 where statistical noise in anyhow very small 
responses may cause those to have negative values. The latter appear to hardly affect the 
extracted parameters. 

Next, the above source of information does not contain numerical data and ER for the 
lower q data . As a consequence we had to limit ourselves to a small data base which is 
bound to influence the FSI parameter functions which increases with decreasing q. 

Our results for T=2.3K are assembled in Table I. We entered in column 2 theoretical 
values for the parameters, calculated as indicated , or set to when impossible to evaluate 
reliably. Notice that the negatives, obtained from a limited s— range, will generate a non 
bound R(q,s), barring a convergent Fourier transform of R(q,s). Actually we did not 
consider ze in the fits and restrict ourselves to a maximum of seven parameters as is discussed 
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below. 

Column 3 is the result of a 7-parameter fit for cumulant coefficients functions with the- 
oretical g-dependence and encompasses therefore all ten data sets in our g-range. Columns 
4 and 5 are 6-, respectively 5-parameter fits like in column 3, when first q*f3^{q) alone and 
then also q* 2 Pi(q) are held at their theoretical values. 

We now reach the extraction of j2 n {q) for each q and the determination, as opposed 
to the above assumption, of their g-dependence. This appeared to be impossible for the 
limited data set available to us. However, Glyde et al had access to far more data and we 
enter their results in the last column taken from Table I of Ref.0. As the reported Z4 is 
zero, a positive value of z 6 is required to give a convergent Fourier transform of R(q, s). 

We hinted in the Introduction, that the misfit of FSI parameters may reflect on the 
meaning of the parameters from which one builds the momentum distribution n(p). In fact, 
the above was the main reason why we repeated the fit of the cumulant expansion coefficients. 
Table I shows even smaller values for z% = q*(3z(q) than reported in00. Although the fitted 
Z3 appears strongly correlated to other parameters, we remained worried by the misfit. The 
above observation led us to ponder about the very possibility to obtain reliable values for 
the cumulant coefficients for the response. 

In fact, we found the results of column 5 in Table I most telling. Fixing the dominant 
FSI parameters z 3t4: at their theoretical values, one expects a 4 , a 6 to settle close to their 
starting values. This appears not to be the case. In fact, the above happens to produce the 
worst agreement between the reconstructed and computed n(p) to be presented below. 

To understand the above, we turn again to the calculated Cl(q, s) which we recall, pro- 
duces very good fits to the data. Although not called for in calculations using a dynamic 
theory, a comparison of ( |5b| ) and ( |20c| ) shows that the expansion of Cl(q, s) in s produces in 



principle the cumulant coefficients and in particular the FSI parameters (3 n {q). Examples 
are O and (pig 1 ). 

The complete cumulant series is of course equivalent to the exact fl, but a truncation 
at some n obviously reproduces behavior up to some relatively low s. The crucial question 
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is to what order one should go and the answer clearly depend on the effective s range of 
each of the component factors. Those are according to (|5b| ), (|IBD and ( |2T|) the single-atom 
density distribution pi(0, s)/p, the FSI interaction factor R(q, s) (or Q(q, s)) and the Fourier 
transform E(q, s) of the ER function. 

In the tests to be reported we compare not data, but theoretical values of the input factors 
pi(0, s), R, Cl and their cumulant approximations R cu,n , Cl cu,n to order n. Fig. 5a shows over 
the relevant s-range reasonable adequacy of the cumulant expansion for pi(0,s)/p to order 
n = 6. Next, Figs. 5b, c give Re and Im parts of fi(g, s), respectively R(q,s) q = 23A 
and shows that the FSI cumulant series, truncated at n — 5, rapidly falls short from of 
the computed functions for s > 1A . Inclusion of the above mentioned, well- determined 
zq ps ;z 64 affects only Ref2, (cf. Eq. ( [20b|) . It extends the agreement between the calculated 



and cumulant expansions for the two FSI functions over a modest additional s-range, but 
does not prevent rapid deterioration of R(q, s), &(q, s) from 1.2 A 1 on. As mentioned before 
the estimated z$ gives rise to a non-bound R(q, s). Since terms in the cumulant expansions 
(p^) have alternating signs of powers in s, it requires knowledge of higher order (3 n in order 
to establish the proper order of the expansion for given s range. 

We return to the interpretation of our results. The cited parameters for finite— n expan- 
sions indeed reproduce very good fits to data, but have no resemblance to the values of the 
parameters from actual (in case theoretical) FSI factors R(q, s), fl(q, s). The failure is due to 
the insufficiency of low order cumulant representations R cu,n , Q cu,n and will only disappear 
if for sufficiently high n, FSI functions will coincide with their cumulant expansion out to 
the relevant medium s. 

The above is to our opinion the source of the discrepancy between calculated and ex- 
tracted FSI parameters. In spite of correlations between fitted parameters, there is sufficient 
meaning to the extracted FSI to justify the above conclusion. Were it not for the overwhelm- 
ing role of the asymptotic part and its representation by a truncated cumulant series, one 
could not trust the extracted a n . 

The sensitivity of the fitting procedure can be judged from Fig. 6, where we reconstruct 
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single-atom momentum distributions from the various a n . Dense dots, short dashes and 
spaced dots are results using columns 3-5 from Table I. Long dashes are from the last 
column, i.e. the fit of Glyde et al. All fits are roughly of comparable quality when judged 
as to their capacity to reproduce the computed distribution n(p,T = 2.5K)0, marked by 
circles. We recall that the theoretical n(p) has been used to compute a n in column 2 which 
produces the drawn curve in Fig. 6 . The moderate misfit for large s of the finite order 
cumulant expansion for pi(0, s) explains deviations from the underlying n(p). 



V. SUMMARY AND CONCLUSIONS. 

We addressed above three topics regarding the response of liquid 4 He, retrieved from 
inclusive scattering of neutrons. Using dynamics we first made predictions for the recent 
MARI data, taken at temperatures both below and above T c . We analysed the range of 
momentum transfers 21 < g(A ) < 29, for which we had available data and experimental 
resolution and obtained good agreement with experiment. 

A second topic was the approach of the response to its asymptotic limit in q for fixed 
scaling variable y. For smooth interactions between constituents, that approach is rigorously 
oc 1/q, but the same is not guaranteed when a strong short-range repulsion is present in the 
atom-atom interaction V. We investigated theoretically the response for q < 300A for 
the actual V with its short range repulsion and found re-confirmed our findings from a few 
years ago based on medium-g data: Final State Interaction contributions, over and above 
the asymptotic part still decrease approximately as 1/q. 

Regarding the g-values for which the response has been measured we repeat what was 
already evident from the older data!: there is no additional information to be retrieved by 
increasing q by less than 20-30% and no new information at all at very high q. 

Our last topic was a re-fit of the expansion coefficient functions, which parametrize data. 
A previous analysis led to a single atom momentum distribution in good agreement with 
computed n(p,T), but did not produce the main FSI coefficient function. Its influence in 
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the large-g region is too marginal, to be extracted from the data with present accuracies. 
This does not change our judgement that little can be added to understanding the data on 
the response of liquid 4 He at high momentum. 

In spite of minor discrepancies, a summary of the treatment of the first two topics is 
definitely positive. Presumably for no atomic, molecular, nuclear or sub-nuclear system 
for which the response has been measured, has one reached asymptotia as clearly and as 
well understood as for liquid 4 He. Of course 'asymptotia' is not simply the mathematical 
limit q — > oo for fixed y. Increasing q requires increasing beam energy, ultimately beyond 
the ionization energy w 39 eV. A description of the response then requires the inclusion of 
additional electronic degrees of freedom to translational ones. 

Finally we venture an outlook for future explorations of inclusive scattering which almost 
certainly implies extension of experiments to larger momentum transfers q. Our judgement 
is based on Eqs. (|5b]) and (|7|). In that representation for the response, the FSI factor R(q, s) 
does not depend on the single particle density pi(0, s) and in particular not on the condensate 
fraction for T < T c . Together with OTP] ) this implies the universality of the asymptotic limit 
(ib~) under all circumstances. One then tends to believe that present experiments and theory 
seem to have exhausted the search for information, contained in the response of liquid 4 He 
at high momentum. 
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Figure Captions 

Figs. la-d. Predictions for the response (|5|), (0) of liquid 4 He at T = 2.3K for q = 
21, 23, 25, 29 A" from (g) using the ladder approximation (|H]) for binary collisions. Those 
results have been folded in experimental resolutions from Azuah's thesis. The same is the 
source of dat J. 

Figs 2a-d. Same as Figs, la-d, for T = 1.6K, computed for rio(T = 1.6K)=0.087. 

Fig. 3a. (fi even (q, y), part of the response even in y for 21 < q(A ) < 300. For q > 25A 
those can not be distinguished from the asymptotic limit F (y), Eq. 

Fig. 3b. q(p odd (q,y), q times the response, odd in y for 21 < q(A 1 ) < 300. 

Fig. 4 . The FSI function R(q = 23A~\y) in fl5aj), computed with the GRS and HNC 
choices for p2, C2 as well as for other descriptions. The drawn line is our result. Long dashes, 
short dashes and dots are from referencesllBi 

Figs. 5. a) Comparison of of the single atom density matrix pi(0, s)/p and its 6th order 
cumulant expansion, b) Real and imaginary parts of the FSI phase s) (q = 2lA 1 ) 
from theory and represented by a 5-th, respectively 6-th order cumulant expansion; c) the 
same for the FSI factor R(q, s). 

Fig. 6. The single atom momentum distribution for T=2.3 K, reconstructed from the 
various cumulant fits, assembled in Table I, including the one from theoretical starting values. 
The curves correspond for increasing column number of the column in Table I including the 
fit by Glyde et al. The circles are values calculated in@. 
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Table I 



Cumulant 


computed 


7-parameter 


6-parameter 


5-parameter 


Glyde et al. 


coefficient 


starting 


fit for prescribed 


fit; fixed 


fit; fixed 


n FT71 




values 


g-dependence 




q*(3 3l q* 2 (3 4 




a 2 (A ) 


0.916 


0.910 


0.913 


0.914 


0.897 


a 4 (A ) 


0.470 


0.553 


0.594 


0.781 


0.46 


a 6 (A ) 


0.337 


0.535 


0.613 


0.700 


0.38 


9*&(A ) 


0.555 


0.237 


0.555 


0.555 


0.33 


g* 2 &(A" 4 ) 


-2.268 


-0.698 


-0.993 


-2.268 





g^54(A" 5 ) 





0.416 


0.851 


0.615 





g* 3 & 2 (A" 6 ) 





-0.152 


1.32 


3.23 


0.201 


g* 2 ^64(A" 6 ) 


» -31.0 








1.539 



Cumulant coefficient functions from data at T=2.3K. The second column gives theoretical 
values. The third column gives 7-parameter fits with prescribed q* behavior from (renor- 
malized) data. Column 4 give fits when q*j3%{q) is fixed at its starting values. Column 5 are 
fits if in addition q^fi^q) is fixed. In the last column are results by Glyde et al, who fitted 
Hn(q) from a large set of data for different q. 
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